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One of the most dramatic results from the first RHIC run are the STAR 
, results for 7r ± 7r ± interferometry. They showed that the ratio of the so-called 

Rout and R s id radii seem to decrease below unity for increasing transverse 
momentum of the pair (Kt)- This was subsequently confirmed by PHENIX 
^sO • which also extended the Kt range of the measurements. We consider here 

the effects of opacity of the nuclei on this ratio, and find that such a small 
value is consistent with surface emission from an opaque source. 

o 
in 

^ : I. INTRODUCTION 

p 

The striking result on tt^tt^ interferometry data obtained by STAR Collaboration |1[ on 
the ratio R out /R S id is not yet understood. This is the ratio of the outward radius, measured 
in the correlation along the direction of the sum of the particles' momenta, divided by the 
sideward radius, reflecting the correlation perpendicular to the beam and perpendicular 
to the outward direction. This result was subsequently confirmed by PHENIX M. In an 
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attempt to understand this result we consider in detail a generalization of the model used 
by Heiselberg and Vischer || to include effects of source opacity at CERN energies, later 
followed by Heinz and Tomasik Q. 

The origin of the problem can be understood from the simplest model of HBT |5|,|| . In 
this model we let particles be emitted from the entire volume of the system. In this case, 
the spatial size which is probed by R out and R S id is of the same order of magnitude, the 
radius of the system. However, R out also has a spatial correlation built in on account of the 
time differences of different emissions 0. R out is the size scale measured along the particle's 
direction of motion. Particles emitted at much different times end up spatially separated 
from one another in this direction. We therefore expect that since the time separations are 
of the order of the size of the system (which for practical purposes seems to be a reasonable 
approximation for order of magnitude estimates at RHIC), then R out / R si( i > 1. Most simple 
models incorporating the above physics give the ratio to be about R ou t/Rsid ~ 2. 

The assumption that the particles are emitted from the entire volume of the system is 
probably a bad one. Data on distributions of particles in STAR and PHENIX || 
suggest that the system is opaque to particles up to large transverse momenta. A more 
reasonable description might be blackbody surface emission. 
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To implement opacity, we consider a model where the matter is emitting from a surface 
at a fixed radius. The system is allowed to undergo 1 + 1 dimensional longitudinal expansion. 
The details of this model and its results are described in the fourth section. Many of the 
features of the model we propose are embodied in the hydrodynamic computations of Heinz 
and Kolb. || The essential difference lies in the treatment of surface emission. 

An important ingredient of this model is decoupling, which we assume occurs at a well 
defined temperature. In the surface emission model we use, this requires that, at the end of 
the process, the decoupling to occur throughout the transverse volume of the system at a 
well defined time. It turns out that, although it is a relevant ingredient (i.e., it contributes 
to the averaging over different emission radii), this is not so important in our computations 
because the surface emission dominates. 

In the second section, before turning to explicit model computations, we review the 
method of the Covariant Current Ensemble formalism for computing two particle HBT 
correlation functions. In this section, we try to be general enough to include the physics 
needed for our various computations. 

In the third section, we consider a simpler version of the above model. We compute the 
various radii for the case of a cylinder emitting at a constant temperature T for a finite 
time At, without introducing any radial flow. We will consider both the case of an opaque 
source and a transparent one. In this oversimplified model, we show that in the opaque 
source limit, that R out ~ J 'At 2 (Kx / ' Ek) 2 + (0.2_Rt) 2 and R S id ~ Rt- For the realistic case 
where At ~ R T , R out /R sid ~ At(K T /E K )/R T , and this ratio can be less than one. In the 

transparent limit R out = J 'At 2 (i£r/ Ek) 2 + R\ and R S id ~ Rt, this ratio is always larger 
than 1. 

In the heavy ion experiments, the ratio R ou t/R S id is larger than 1 at AGS to SPS energies 
for 7r _ 7r~ pairs ||||. For it is also larger than 1 at AGS (R ou t/Rsid ~ 1-3 - E859 Collab. 

Hill)) but slightly smaller than 1 at SPS, according to NA44 data f§. 

Although the value of R out /R S id is above one at the AGS, and it is relatively close to one 
at the SPS, 0- |L(J, the data at RHIC gives an even smaller value. This suggests that at 
RHIC this ratio reflects the opacity of the emitting surface. The origin of the smaller value of 
Rout/Rsid for the opaque case is due to two effects. First, the system radiates from a smaller 
geometrical region than in the transparent case. Second, the emission itself is preferentially 
from the region of the cylinder closest to the radial vector made by R ou t, where this vector 
is a radial vector in the direction of the pair transverse momentum K^. This is because 
radiation from this region is normal to the surface, and there is a factor of cos0 associated 
with the flux in the direction of R ou t, where <fi is the angle between H ou t and a normal vector 
to the surface. Finally there is also the emission from the decoupling volume which occurs 
at a well defined time, which we shall see does not significantly contribute to the flux of 
particles, although it still contributes to the averaging over different emission radii. 

The case of 1 + 1 dimensional longitudinal expansion is more complicated. In later 
sections, we shall use a simple model of the matter, where there is a real first order phase 
transition between a parton gas (gluons and massless quarks) and a pion gas. The energy 
is emitted from the partonic phase. This quark and gluon matter is assumed to be directly 
converted into a flux of pions with the same energy and a blackbody distribution at the 
temperature of emission. This energy conservation condition allows us to directly take a flux 
of gluons and quarks and convert it into a spectrum of pions. It of course will mess up details 
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of the fragmentation, and generates some increase in entropy, but for our purposes such a 
crude treatment is sufficient to demonstrate the physics of opacity in a semi-quantitative 
way. The emission of the hadron degrees of freedom is taken to be entirely pions. We will 
find that the dominant emission at RHIC energies comes from times after the beginning of 
the mixed phase, but there is a significant mixture of gluonic and quark radiation which 
must convert into pions. 

The characteristic time scale for emission is of the order of the radius of the system 
before the radiation is significantly attenuated. This time scale is of the order of the size 
of the system for RHIC energies. Our model which incorporates this physics is described in 
detail in Sec. 4. In the fifth section, we describe some of the details of the computation of 
the HBT correlation function. 

In the final section, we discuss the limitations of our approach. The severest is that if 
we try to predict the values of HBT radii as a function of Kx of the pair, we do not get the 
correct transverse momentum dependence. We do however get the correct dependence 
for R out / R S id- 



II. THE CO VARIANT CURRENT ENSEMBLE INPUTS 

To compute the emitted spectrum and the two particle distribution function, leading to 
the interferometry relations in which we are interested, we consider the Covariant Current 



Ensemble formalism, [14], fT5|| . In this formalism, the two particle correlation function can 
be written as 



P2(h,k 2 ) , \G(h,k 



° {ku h) P x {k\)P x {k 2 ) 1 + G(k h h)G(h, k 2 ) (1) 

where P\{ki) and P 2 (ki, k 2 ) are, respectively, the single particle distribution and the proba- 
bility for simultaneous observation of two particles with momenta k\ and k 2 . The average 
momentum of the pair is defined as K = (k x + k 2 )/2, and the relative momentum, as 
q = ki - k 2 . 

The complex amplitude, G(ki,k 2 ), is written as 

G(h, k 2 ) = [d 4 p [ d 4 x eV" D(x,p) j*(^£)j (^P) , (2) 
J J mm 

and the single-inclusive distribution, G(k i} ki), is written as 

G(h, ki) = J d4 P^P)\U^)\ 2 = ^^=E l ^ = P 1 (k l ) . (3) 
We see that 

D{q,p)= f d 4 xe iq ^ D{x,p) , (4) 
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where D(x,p) is the normalized phase-space distribution at the instant of the emission. 

One should note that there is some arbitrariness in how we decompose into D(x,p) 
and the densities jo- We shall follow the original proposition in ||14|| , |Tj| , and choose to 
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include the thermal aspects of the collision in jo, and the geometrical plus collective effects 
in D{x,p). In any case both will be defined for the problems at hand. 
There are several different possibilities for D(x,p): 

The first is radiation from a transparent cylinder which is not expanding in the transverse 
directions and has no longitudinal expansion. 

D(x,p) = 9v k exp(-t 2 /2At 2 ) 5(E p - p°) 5 2 (p T )5(p z ) . (5) 

In this expression p is to be interpreted as the momentum of a particle at rest in the 
matter from which it is emitted, and reflects the collective motion of the system. For matter 
at rest, p = 0. The factor g p counts the number of degrees of freedom of particle being 
emitted. The factor of e~ l / 2A * is a Gaussian parameterization of the source emission rate, 
so that the characteristic emission time is At. (The blackbody rate emerges naturally in 
this picture, and k is a factor which modifies beyond the blackbody rate). 

The second case is emission from the surface of a cylinder without longitudinal expansion. 
Here we must account for the fact that the particle is emitted with a flux factor which 
depends upon the angle between the normal vector to the surface and the direction of the 
pair transverse momentum K?- There is also a factor of <5(ty — Rt) which requires that 
particles are emitted from the surface and a factor of G(cos0) which requires they come 
from the side of the surface from which they are emitted. 

D(x,p) = g p n exp(-t 2 /2At 2 ) 5{r T - R T ) S(E p -p°) <5 2 (p T ) 6(p z ) cos0 9(cos0) . (6) 

Finally, there is the case of emission from the surface of a cylinder with longitudinal 
expansion 

D(x,p) =g p n5(r T -R T ) 5(E p -p°) 5 2 (p T ) <% - r?) cos0 6(cos0) , (7) 
The factor of 5(y — if) is the correlation between the velocity and coordinate of the emitting 



surface built into the Bjorken model [|12[]. The factor of k controls the rate of emission from 
the surface as a function of time. We will fix this by requiring the system to be either a 
Quark Gluon Plasma, a hadron gas or a mixture, with blackbody radiation and the mixture 
determined by thermodynamics. 

The currents correspond to thermal distributions and are written as 

Here = p^/ m is the four vector flow velocity which is the four velocity of the emitting 
surface. In the case where the currents are connected to the models illustrated by Eq. (H)-© 
above, u M = (1,0), i.e., associated to the absence of expansion along the longitudinal and 
transverse directions in those cases. We indicate by T(r) the dependence of the temperature 



in the proper time, r, prior to the beginning of the mixed phase, as given by Eq.fll7l), in 
Section IV. 

In the ideal Bjorken picture with no transverse flow, the 4- velocity and the momentum 
of the emitted particle, k^, can be written as 

= (cosh?7, 0, 0, sinh?]) ; = (uit coshy, k T , ttit sinhy) . (9) 
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III. A VERY SIMPLE MODEL FOR R ou t/Rsid 

We first work out the simplest model for emission so that we can get some conceptual 
understanding of the physics involved in the various cases. We shall consider a system 
without radial flow. The approximations will be fixed in the next section when we consider 
a more realistic model which we compare with data. 

First consider emission from a transparent cylinder. A little algebra gives 



Cout = 1 + exp 



(K 2 + m 2 ) 



[2Ji(qoutRT)/(q utR T )f 



(10) 



and 



a id = 1 + [2Ji(g S idiM/teid#T)f 



(11) 



where J±(x) is the Bessel function of the first kind and order one. 

In this equation, it is absolutely clear that R out > R S id, and that if At ~ R, then 

Rout ~ 2R.sid 

Now let us consider the result for surface emission. In this little algebra yields 



where 



On the other hand, 



Cout = 1 + exp 



1 f n 

(x) — — / d0sin(0) expfecsin 

2 Jo 



(K 2 + m 2 ) 



IiqoutRr) 



a 



sid 



1 + 



sin(q sid R T ) 
(qsidRr) 



(12) 



(13) 



(14) 



In Fig. 1, we plot C si d and, for stressing the geometrical differences only, we ignored the 
time contribution to C out , by fixing At = in the plots. Fitting the curves by Gaussian 
distributions, we find that, in the opaque case, R e JJ ~ Rt x 0.53 and R e J u { = Rt x 0.22. In 
the transparent case, the fit results in an effective radius of R^{ ~ 0.61-Rt- 

The basic result we get from this analysis is that for the opaque cylinder, unlike the 
transparent one, we can easily have R out < R S id- In fact for R out ~ At(K T /E K ) (since 
the R T contribution is very small), the ratio of R ou t/R S id ~ At(K T / E K ) / R T unlike the 
transparent cylinder case where R ut/R Si d ~ \j At 2 (K T / E K ) 2 + R 2 T /R T . 



IV. A SIMPLE MODEL FOR R ou t/R s id 

In this model, we incorporate longitudinal expansion. Here the system cools as it ex- 
pands, so we need to have the rate of emission vary as a function of time. 

We need to have dynamical description of the microphysics to be able to do this. Our 
more or less conventional model is a Quark-Gluon Plasma (QGP) phase at some temperature 
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above the critical temperature, T > T c . We choose this critical temperature T c = 175 MeV 
to be consistent with lattice Monte-Carlo data. The system goes into a mixed phase of QGP 
and hadron gas at the critical temperature, and we take this hadronic gas to be composed 
of an ideal gas of pions. Below T c , it is only an ideal gas of pions. While in the low 
temperature pionic phase, the system cools until it reaches a decoupling temperature, Tf. 
We choose Tf = 150 MeV to be consistent with the typical energy per particle observed in 
the RHIC experiments, and to fit the observed pt distributions of pions. However, these 
values chosen for T c and Tf are not crucial for qualitatively reproducing the results we discuss 
here, since these results are only weakly sensitive to these particular values. 

The system produced in a heavy ion collision will expand and the temperature will 
gradually decrease. The initial expansion is in the Quark-Gluon Plasma phase, and the 
system expands longitudinally. After this initial stage, lasting (r c — tq), the transition 
temperature, T c , is reached and the evolution continues in the mixed phase, during which 
the temperature remains constant with time. The mixed phase continues for a longer period, 
ending after an elapsed interval {jh — r c ). Then, the system converted into a gas of pions, 
further expands until the decoupling temperature, Tf, is reached. At this point, the system 
is quite dilute, and much of the particles have been evaporated from its surface. Thus, 
we consider that, once Tf is reached, the system decouples in an instantaneous volumetric 
emission. 



The Bjorken hydrodynamical model [12| should be able to describe the system during 
its evolution from formation until the time it breaks up. We will supplement this with 
radiation from the surface of the matter. We will take the radius at which this radiation 
takes place to be a constant and equal to the radius of the nuclei. We only consider impact 
parameter zero collisions. In the last stage of the system evolution we consider a volumetric 
emission at freeze-out but, as in the surface emission, no transverse flow was introduced in 
the computation. This could be one of the main reasons for obtaining a much weaker Kt 
decrease of the transverse radii, as compared to RHIC data [|T],0. 

Another ingredient in our model is the hypothesis that the system will emit from its 
external surface similarly to a blackbody, starting shortly after being formed, at r = r . In 
this way, quark and gluon degrees of freedom have to be considered in the QGP and mixed 
phases. The hadronic degrees of freedom should, in principle, include a complete set of 
resonances later decaying into pions. However, in this initial description, and for the sake of 
simplicity, we will consider a hadronic gas constituted of only pions. No complex mechanism 
for the QGP hadronization will be considered in detail at this point, although hadronization 
must take place. In other word, in first approximation, we will consider the evaporation of 
"gluons" and "quarks" (as hadronized pions) from the external surface of the system in the 
same way as emission of pions, except for the number of degrees of freedom. 

We can estimate the emitted energy as well as the total entropy associated to each stage. 
In the initial phase, lasting from r to r c , we can estimate the emitted energy as a function of 
time considering the emission by an expanding cylinder of transverse radius Rt and length 
h, in a certain time interval between r and r + dr by 

dE in = —koT^ 2vr R T h dr - - cxT 4 tt R 2 T dh , (15) 

3 

where the first term comes from the blackbody type of energy radiated from the surface of 
the cylinder, and the second term results from the mechanical work due to its expansion. 
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The k factor was introduced to take into account that the system has some opacity to 
surface emission. The constant a is is proportional to the number of degrees of freedom in 
the system. 

By integrating Eq. (^) we get for the energy density (i.e., e = E/V) 

tin = eo(— ) 3 e r t k . (16) 



From the above expression we see that we obtain an extra multiplicative factor, 
e r t {t T0> , in addition to that coming from the Bjorken picture. Remembering that, in 
the Bjorken picture, the relation between the energy density, e, and the proper-time, r, 
is given by e/eo = (ro/r) 4 / 3 , and that e = crT 4 for a blackbody-type radiation, then the 
variation of the temperature in the initial stage, i.e., prior to the beginning of the phase 
transition, follows immediately as 

T(r) =T (^ e ^ (T - ro) . (17) 

T 

In order to fix the initial conditions of the evaporating and expanding fireball, we follow 
the evolution of the entropy, using the observed final particle multiplicity density dN/dy 
and its relation to the final entropy density dS/dy as a constraint. We decompose the total 
entropy S to t into its contributions Si n from the fireball interior and from the emitted particles 
Semit, all of which are functions of the proper-time. During the first stage, from tq to the 
beginning of the mixed phase at time r c we, can write 

= % = ke + p)* = so(-)e-^ (T - ro) - S m = S e-^ (T - T0) , (18) 

VI T 

where s = f (o"e 3 ) ) 1/ ' 4 is the initial entropy density. In the second step of Eq. fll8|) we used 
that, for the Bjorken type of longitudinal expansion, the volume V(r) increases linearly with 
time. 

On the other hand, the entropy associated with the emission can also be estimated. For 
this, we write the element of emitted entropy, dS em it, in the interval between r and r + dr 
as proportional to the entropy on the surface (i.e., to the entropy density times the surface 
area, 2tiR t t). The emissivity, k, is the proportionality constant. If we divide both sides by 
the corresponding volume, V = tiR^t, we then get 

^ = ^s m dr = ^s (^)e-^-^dr - S emU = ts [l - e"*^] . (19) 
V ii r p JrC r r T 6 

This results in the total entropy of the initial stage as being 

S tot = p - ifib e-^ {T - To) , (20) 

so that, at r = r — » Stot = So- Note that this result requires that there be entropy produced 
during the emission from the surface with the added input that the temperature changes as 
a function of time. In the mixed phase there will be no such effect, since there the emitted 
quanta are simply the number of the quanta in the system, and no entropy increase in particle 
number is generated. Note that the entropy for a massless gas is directly proportional to 
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particle number. The situation is simply different when one has an expanding system with 
a variable temperature. 

During the phase transition, the energy and entropy can be estimated similarly, leading 

to 

E emit = 4(1 - e-^ (T - Tc) ) ; S m = S c e -^ {T - Tc) ; S emit = S C {1 - e"^^) . (21) 



These results follow simply from Eq. ( 15|) when the second term on the right hand side is 
set to zero. During the mixed phase no work is done in expansion because the temperature 
does not change and the expansion conserves entropy. Therefore only the first term on the 
right hand side of Eq. ( |I5|) contributes, and it is easy to integrate, since it is 2kE(1t/R. 

From the above relations we can see that, during the phase transition, the total entropy, 
Stot = S in + Semit = S c is conserved. During this extended period, the temperature remains 
constant with time (T = T c ), so that the previous relation in Eq. QTTj ) no longer holds. Note 
that there is a difference in the exponential behavior in the mixed phase relative to that in 
the QGP phase. This is a consequence of the fact that the temperature is time dependent 
in the QGP phase, but time independent in the mixed phase, so that in the first case the 
rate of emission from the surface is different. 

During the phase transition the system is in a mixed phase of Quark-Gluon Plasma and 

hadronic gas. If the fraction of the fluid in the QGP phase is /, then in the hadronic (pion) 

phase it would be (1 — /). On the other hand, from Eq. (pip, we see that the portion of the 

2k i 



entropy density still in the system is given by s in = s c (— ) exp [— jt-(t — r c )]. Consequently, 



Sin = f s QG p + (1 - /) s h , where s in (r c ) = s c and s(r h ) = s c (^) exp [-^(r h - r c )] . If we 
substitute these expressions into that for s in , we then get 

f= { r h -r c J r ' (1 " /} = \ Th _ Tc e -*(— > ) r ' (22) 

Finally, we need to estimate the initial values T and r , as well as the proper time, 
r c , corresponding to the on-set of the phase transition. We estimate r by means of the 
Uncertainty Principle, i.e., {E )tq m h, and by 

{Eo) ~ Jdpp* e-v/T, ~ 3To ' ( 23 ) 
from which we easily determine the initial time as 

h °- 197 f / ,nA\ 

T °~3T = 3(T /GeV) fm/C " (24) 

On the other hand, the constraint on To and tq has to come from the experiment. At 
RHIC, the average produced pion multiplicity per unit of rapidity is M ~ 1000, which should 
be proportional to the initial entropy, So, i.e., 



S = TAf 



G?g + S q )x(^T 3 



ttR^tq . (25) 



In the above expression, we have related the initial entropy, So = SoVo to the initial 
entropy density and the initial volume of the system, Vq = ttB^To, this last one estimated in 
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the Bjorken fashion. The degeneracy factors, g, are given by the gluon degrees of freedom, 
g g = 2{spin) x 8(color), and the quark/anti-quark degrees of freedom, g q = |[2(spm) x 
2(9 + o) x 3(color) x Nf(flavor)], which add up to g qgp = g g + g q . In the case of pions, the 
degeneracy factor is g w = 3. 

From Eq.(p4|) and (53) we can determine Tq as 



T = VA/T 



Og + 0q) X 



4tt 3 

270 (0.197) 1 



-1/2 

GeV . (26) 



As an example, if we take T = 3.6, as estimated by the entropy per particle (S n /N n ) of a 
pion gas at freeze-out, then To ~ 411 MeV and tq ~ 0.160 fm. 

For estimating the instant corresponding to the beginning of the mixed phase, r c , we 
consider Eq.([T7|) at r = r c , resulting in 



3k ^ / 1 1\ \ ^ 3k „_ 



r c e ^ Tc = (^J r eW T0 , (27) 

which can be numerically estimated for fixed values of k and T c (for this, we will consider 
T c = 175 MeV). 

In order to estimate the instant corresponding to the end of the mixed phase, t/j, we need 
to know s c and §h- For the first one, we consider a system of gluons and (massless) quarks 
forming an ideal gas, resulting in s c = (g g + g q )^-Tj?. For estimating 7^, we consider the 
pions as massless particles, while in the system, leading to §h = g^^-T^. Then, equating 
the expression for s in at r = r h , we obtain 

r h e^ Th = {^j^j t c , (28) 

which can be estimated numerically for a fixed value of k, since r c was already determined 
by Eq. fl27|). 

We assume that, at the end of the phase transition, corresponding to r = r h , the system 
is an ideal gas of pions (no resonances are considered in this initial estimate), which continues 
to expand and cool down, until the temperature drops to T c = 150 MeV, corresponding to 
an instant r = r/. At this point, whatever is remnant of the system breaks up. The portion 
of the pions still in the system at that time we call V, i.e., the fraction of the system that 
is emitted from the entire volume at the time r/. This decoupling instant can be estimated 
by an expression similar to Eq . (fl7|) ) , with T replaced by T c , and r by r c , resulting in 

(29) 

which can be numerically estimated for fixed values of k, T c , and Tf, and corresponding t^. 

We illustrate in Table 1 below these variables for two different assumptions on the emis- 
sivity, k. The values of the temperatures considered were To ~ 411 MeV (obtained from 
Eq. (Ej|)), T c = 175 MeV and Tf = 150 MeV. We also include the corresponding estimates 
of the fraction of particles emitted from the surface, S, in the interval To < r < Tf, as well 
as the fraction from the volumetric emission, V, at r = Tf. 
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To estimate the fraction of the particles emitted from the surface, S, and from the 
volume, V, we proceed as follows, estimating the contribution from each stage of the system 
evolution. For that, we estimate the emitted entropy in each stage, similarly to the procedure 
described in Eq.flTD-(|TD. 

The fraction of the input particles, Af, emitted in the first period, r < r < r c , is given 

by 

$ = *(l_ e -*C~->). (30) 

Similarly, the fraction A/^/A/" emitted during the phase transition, r c < r < th, can be 
written as 

Af 2 _3^6, ~\ la 

Af 



qi = e -afcfc.-*0(i _ e -#^-^)) . (31) 



The fraction emitted during the pure pion phase, < r < Tf, up to reaching Tf, can be 
estimated by 

^£ = 1 e -27^( T <=- r o) e~^F (Th ~ rc) (l - p-m^^f-^h f32) 
AT 3 v ; v ; 

Then, the fraction of the particles emitted from the surface up to r = Tf is 

= • (33) 



Finally, the remnant fraction at r = Tf, emitted from the entire volume, is given by 

(34) 



V _ ^ _ -^(r c -r ) -|^(r fc -r c ) -^( r/ -r fe ) 



Af Af 

As a result, the ratio of the total number of emitted particles, Atot, with respect to the 
input number, Af, is given by 



A/tot l_ 
AT Af 



(M+A^+ATs+A^) • (35) 



TABLE 1: Values of the proper-time parameters r , r c , r^, 77, as well as of the 
surface, S (for r < r < ry), and the volume, V (at Tf), emitted fluxes, for two 

values of the emissivity, k. 





TO 


T c 


Th 


T f 


S/N 


V/Af 




(fm/c) 


(fm/c) 


(fm/c) 


(fm/c) 


(t < r < r/) 


(at r/) 


1 


0.160 


1.54 


5.73 


6.97 


0.844 


0.156 


0.5 


0.160 


1.75 


8.37 


10.5 


0.758 


0.242 
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As we can see from the fractions in Table 1 above, in our model there is a small increment 
(~ 7.6 — 10%) in the total number of particles with respect to the initial one. We show in Fig. 
2 the evolution of the emitted flux with proper time, by plotting the fractions normalized to 
the total number of produced pions, as a function of r. The two different cases correspond 
to different surface emissivities. The curves end at the decoupling temperature. In the case 
of the greatest rate of surface emission, about 84% of the total radiation comes from the 
surface and the rest from the decoupling volume. It is about 76% for the other case we show 
for illustration. In the first case, the system decouples at a time of about 7 fm/c and in 
the second case at about 10.5 fm/c. We assume a decoupling temperature of 150 MeV. We 
should note that the decoupling proper-times within our model are significantly shorter than 
those in hydrodynamics. This is mainly due to the fact that we allowed a higher surface 
emissivity than those kind of models, which may be the key to explain the proper-times 
observed at RHIC, which are smaller than expected from hydro predictions. 



V. SINGLE- AND TWO-PARTICLE PROBABILITY DISTRIBUTIONS 

We define the average momentum of the pair as K = h(ki + k 2 ) and the relative momen- 
tum as q = (ki — k 2 ). They satisfy q^K^ = 0, and, consequently, the temporal component 
of q 11 can be written as q° = In the limit that is interesting for interferometry, we can 

consider |q| <C |K|, which implies that K° « y^K] 2 + m? = Eg. 

In both expressions for C(ki,k2) and C(ki,ki), we see that, due to the form of the 
phase-space distribution in Eq.(^), integration over the variables involving delta functions 
are straightforward. We should remember that, due to the factor O(cos0) in Eq.flT]), the 
integration over dcj) runs in the interval [— 7r/2, tt/2], while the limits on the proper time 
integration would be [tq,tj]. However, as the system has different composition in each 
phase, we should split this time integration to be [t ,t c ], [T c ,r h ], and then [r^Ty]. The 
rapidity integration should run, in the Bjorken picture, from (—00, +00). 

We recall that / and (1 — /) are, respectively, the fractions of the system in the QGP 
phase and in the hadronic phase, according to the expression given in Eq. (|22|) . Then, by 
taking into account the above observations, the expression for the complex amplitude, after 
some algebraic manipulation, can finally be written as 




dy[Ex cosh y — Kl sinh y] 




exp \-^(E K coshy - K L smhy)(-) 1/3 e K{T - To)/{2RT) } + 




igp 




ir I"eJ It cos a cosh y + t L ( "g^ cosh S/- sinh V)\ -7fr (E K cosh y-K L sinh y) 




, iT lE^lT cos a cosh y+q L ( -gi cos h y-sinh y)] A. (E K cosh y-K L sinh y) 
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f T J j itF-et^-Qt cosacosh«+qr (-5^- cosh y— sinh «)1 

exp [- — cosh y - tf L sinhy)(— ) 1/ V (T_T * )/(2iiT) 

7T 

r % , , K 



} + 



p2tt 




/ ^ 




/o 


Jo J- 



IT 



Jj£ (Jt cos a cosh j/+<?l ( ^ cosh y-sinh y)\ -iq T Rt cos (a-cp) ^~ ^j( e k cosh y-iii, sinh y) (36) 



The five terms composing the correlation function represent, respectively, the emission 
from the quark and gluon initial stage, their contribution during the mixed phase, the pion 
emission also from the surface in that phase, the emission during the pure pionic phase 
up to reaching the freeze-out temperature, Tf, and finally, the volumetric instantaneous 
decoupling once it was reached. 

Similarly, we can write the spectrum, as 
G(ki, hi) oc kRt / <i0cos0 / dy [Ei coshy — ki sinhy] 

J-\ J-oo 

{ <? qg p P t drexpf-i-^ coshj/ - k iL smhy)(-) l ^ T - T ^ 2R ^] + 

g qgp f h TdT (f) e -4-^co^y-k lL sinhy) + ^ T dr ^ _ J) co S hy-fc <£ sinhy) + 

g n [ Tf rdr exp[~(E iC oshy-k iL sinh y )(L)^ e <r-r h )/(2R T)] } + 

J r h ±c T h 

[ dcj) [ r T dr T [ dy [Ei cosh y — k iL sinh y]e T f ^ cosh -y K L svahy ^ 

JO JO J-oo 



7T JO 



We remind that G(ki,ki) is the spectrum as written in Eq.(^|). If we then integrate 
separately the terms of Eq. (|37|) in d 3 ki/Ei, i.e., in the intervals t$ < t < r c , r c < r < th, 
Th < t < Tf, and at r = Tf, we recover the number of emitted particles in each interval, 



as seen in Eq.([3(]), (plf) , (p2|) , and (p4|), respectively, except for an overall normalization 
constant, which is cancelled when we estimate the interferometric relations, as in Eq. ([!]). 

For estimating the (real) amplitudes G(ki,kj) in the denominator of Eq. ([I]), we wrote 
ki = K + q/2 and k2 = K — q/2, from the definition of the momenta K and the relative 
momentum q, since we are not generating the individual momenta and later averaging over 
all of them, as done in the experiment. In this first approach, that is the way we connected 
the momenta hi in the spectra C(k i} ki) to the momenta appearing in the complex amplitude, 
C(k 1 ,k 2 ) = C(q,K). 

This means that we can write 



|k Tli2 | = ^(K T ±q T /2) 2 = JK 2 c + ^-±K T q T cosa . (38) 



\k Ll , 2 \ = \lKl + ^± K L q L ; E l>2 = ^w? + (K T ± q T /2)2 + (K L ± q L /2f . (39) 
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In these equations, a is the angle between Ky and q^- 

Due to the azimuthal symmetry of the problem, we can choose along the x-axis, 
without any loss of generality. In this way, we see that the component of the relative 
momentum in the outward direction (q T || K T ), q TouU will be along the x-axis, while the 
sideward component (q T _L K T ), qT sid , will be directed along the y-axis, i.e., q x = qT out and 
q y = qT aid - This implies that, in the first case, we chose a = and, in the second, a = 7r/2. 

In order to check how the spectra estimated within our model and the above discussed 
relations behave compared to data (PHENIX minimum bias ||13||), we plot the single-inclusive 
distribution in Fig. 3. From here on we limit our estimates and discussions to the central 
rapidity region, i.e., yi — (which implies that k ih = 0, and, consequently, K L = and 
qL — 0). In Fig. 3, we show the spectra, within a constant arbitrary normalization for the set 
of interrelated values shown in the first line of Table 1, and for two values of the emissivity, 
k = 0.5, 1. We see that both curves describe well the spectrum in the low momentum region 
of the pions, up to roughly k iT ~ 1 GeV/c. We shall consider these possibilities also when 
studying the correlation function and fitted radii. 

In Fig. 4 we show the correlation functions C(qx out ) vs. qx out (solid curves), corresponding 
to a = 0, and C(?t sW ) vs. qr std (dashed ones, very close to one another), for a = ir/2, in 
the same plot for better visualize the differences. When calculating the correlation function 
in terms of qT out , we fixed qT sid = (remember that we had already fixed q^ = 0, as a 
simplifying assumption). They are displayed for three values of the average pair momentum, 
Kt = 0.17,0.47, and 0.80 GeV/c, and emissivity k = 1. 

In Fig. 5, we illustrate the behavior of the correlation functions, C(qx out ) vs. qx out (solid 
curves) and C(qT sid ) vs. qr sid (dashed curves), similarly to Fig. 4, but with emissivity k = 0.5. 
Again, the curves were estimated for K T = 0.17,0.47, and 0.80 GeV/c. We see that in this 
case the correlation functions versus qT out are always narrower (consequently, the radii are 
bigger) than those curves versus qr sid - This originates in the higher contribution from the 
volumetric term with respect to the surface ones, coming from their lower emissivity in this 
case, since k = 0.5 in that plot. The larger emission duration is also a consequence of this 
reduced surface emissivity with respect to the previous case discussed in Fig. 4. 

Although not shown for sake of clarity in the plots, the correlation functions, C(qT out ) x 
q Tout and C(qT sid ) x qr sid , were also computed for other three values of K T , i.e., all together, 
for K T =0.17, 0.27, 0.38, 0.47, 0.63, and 0.80 GeV/c. Each of the curves, either for k = 0.5 
or for k = 1, were fitted by Gaussian distributions, in the regions were their behavior 
could be reasonably well approximated to that distribution. In this way, we obtained the 
corresponding average values of R out and R S id, as shown in Fig. 6 and 7, and from those, 
we estimated the ratio R ou t/Rsid- We see that R si d as a function of K T is basically flat 
in both cases, since we have considered no transverse flow in the computations. However, 
our R out decreases with increasing Kt, although not as much as suggested by data. This 
is a consequence of the time dependence of the temperature, that the higher momentum 
particles are emitted at earlier times. The values obtained for the ratio are plotted in Fig. 8 
together with the preliminary STAR (filled triangles) and PHENIX (filled circles) data for 
both 7r + and 7r~. We can see that our results were highly successful in describing both sets 
of data for k = 1, but the curve corresponding to k — 0.5 is away above the data limits, 
suggesting that we should have a high emissivity along the system history in order to explain 
the data trend. 
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VI. SUMMARY AND CONCLUSIONS 



This simple model works well for the ratio of R ou t/R S id and suggests that the origin of 
the experimental value lies in the opacity of source emission and the relatively short time 
of decoupling of the longitudinal expansion, i.e., of the order of the nuclear radius. The 
principal reason why we are able to get such a small ratio of R out /R S ide is probably due 
to a combination of two effects. The first is that the surface is opaque, and whatever is 
emitted from the surface will have a small value of this radius. The second effect is that 
we allow black body radiation by gluons when the surface is very hot. This allows a much 
larger contribution from surface emission than is typical of what happens in hydrodynamical 
simulations where particles are emitted by Cooper-Frye decoupling from a surface at very 
low temperature. In fact we find that about 80% of the emission comes from the surface. 
The fact that so many particle are emitted from the surface at early times also means that 
the longitudinal decoupling time in this computation is significantly shorter than would be 
the case for hydrodynamic simulations, and this again goes in the direction suggested by 
the RHIC data, where the longitudinal time scale did not grow as much as might have been 
expected between SPS and RHIC energies. 

The model also describes the typical source radii reasonably well, but not the Kt de- 
pendence of these radii, as seen in Figs. 6 and 7. This suggests that the time variation 
of the emitting radius and the introduction of transverse flow may play a significant role, 
|Tj| - flT8"| . If there is a time variation of the various radii, this will be correlated with the 
typical momentum scale of emitted particles, since the earlier is the time, the hotter are the 
particles. We are sensitive to such variation since we allow emission from the hot surface at 
early time. Also, a proper treatment of the decoupling is not included in our computations 
and this certainly will affect the results, although it might also suggest modification in the 



treatment of decoupling, ]Tj|, fT"9" |. 

At a minimum, these computations suggest that the problem in describing the various 
HBT ratios lies not so much in R out /R si d as it does in computing the full set of radii and 
obtaining a comprehensive and complete description of all of the above within one dynamical 
model. The ratio of R out /R S id may very well be independent of many of the variations within 
these models since it is dimensionless, where the dimensionful values of the various radii 
are sensitive to changes of time and size scales. Such a complete dynamical computation 
allowing the possibility of surface emission at very early times has not been implemented. 
It may in fact be more complicated than we suggest: Perhaps the surface emissivity is a 
strong function of the momentum of the emitted particles. We certainly expect that high 
momentum particles are more easily emitted than are low momentum ones. 
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FIG. 1. Illustration of the correlation functions in the very simple model for R ou t and R s id as a 
function of the corresponding variable qT out and qr Bid - For emphasizing the geometrical differences, 
we considered At = in the plots of the C(qT out ) vs. qr out - The set of points and the fitting curve 
in the middle correspond to both C(qT out ) vs. qT out and C(qT aid ) vs. qT std hi the transparent case, 
since no time dependence is included in the qT out variable in the above plot. The narrower and 
the wider sets correspond, respectively, to C(qT sid ) vs. qT M and to C{qr out ) vs. qT out m the opaque 
case. We see that, as a result of the opacity of the source, this last set is much broader than the 
first one. 
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FIG. 2. The flux of emitted particles per unit time is shown, normalized to their total number, 
for two values of the emissivity {k = 1 and k = 0.5), at each instant, starting at tq = 0.16 fm/c, 
passing by the beginning of the phase transition, at r = r c , then through its end at r^, and finally, 
stopping at the instant when the pionic system reaches tj . We clearly see the fraction of the pions 
still in the system at that instant, which is then immediately emitted from the entire volume. In 
this example, we fixed T = 411 MeV, T c = 175 MeV, T f = 150 MeV, R T ~ 7 fm/c, and the 
number of quark flavors in the QGP phase to be 2. 
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FIG. 3. The prediction based on our model for the transverse momentum distribution of emitted 
pions is shown. The points are from the minimum-bias data from PHENIX Collaboration. The 
curves correspond to emissivity k = 0.5 and to k = 1, without the inclusion of transverse flow. 
We see that both cases describes data on spectrum well in the low pion momentum region, up to 



about k 
MeV, T c 



« 1 GeV/c. The parameters used are explained in the text, corresponding to To 
175 MeV, Tf = 150 MeV, and the transverse radius, Rt ~ 7 fm/c. 
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FIG. 4. The correlation functions, C(q Tout: K T ) vs. q Tout (solid), and C(q Tsid ,K T ) vs. q TsU 
(dashed), are shown for three distinct values of the average pair momentum, Kt, corresponding to 
the phase transition temperature T c = 175 MeV, Tf = 150 MeV, including a volumetric emission 
when the hadronic (pions) system reaches r = rt. We see that the width of the curves as function 
of qT out increase (or conversely, the radii decrease) with increasing Kt, whereas the curves for 
different qT aid show no visible variation, as would be expected since no transverse flow is considered 
in the computation. The same input parameters were adopted here: To = 411 MeV, T c = 175 
MeV, Tf = 150 MeV, the transverse radius, Rt ~ 7 fm/c, and emissivity k = 1. 
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FIG. 5. The correlation functions, C(q Tout: K T ) vs. q Tout (solid), and C{q Tsid ,K T ) vs. q TsU 
(dashed), are shown for three distinct values of the average pair momentum, Kt, as in Fig. 4, 
but now with reduced emissivity, k = 0.5. We see that the width of the solid curves still increase 
(or conversely, the radii decrease) with increasing Kt, but they are all smaller than the width 
corresponding to the dashed curves. The same input parameters as before were adopted here: 
T = 411 MeV, T c = 175 MeV, T f = 150 MeV, the transverse radius, R T ss 7 fm/c. 
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FIG. 6. Our results for R ou t (outward) radius are shown as a function of the average pair 
momentum, Kt, for the two cases discussed before, corresponding to 50% emissivity and to k = 1, 
without inclusion of transverse flow. The experimental data points from STAR (triangles) and 
PHENIX (circles) are also included in the plot. The values of the parameters are the same as 
in the previous plots, i.e., To = 411 MeV, T c = 175 MeV, Tj = 150 MeV, the transverse radius, 
Rt ~ 7 fm/c. 
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FIG. 7. Analogously to the previous case, we show here the results for R sic i (sideward) radius 
vs. Kt, studied for k = 0.5, 1. No sensitivity to k is seen, since no transverse flow is included 
in the computation. The experimental data points with error bars are from STAR (triangles) 
and PHENIX (circles). The values of the parameters are the same as in the previous plots, i.e., 
T = 411 MeV, T c = 175 MeV, T f = 150 MeV, and the transverse radius, R T » 7 fm/c. 
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FIG. 8. The result corresponding to the ratio Rmit/Rsid of the outward radius by the sideward 
one is shown within our model. We see that the ratio corresponding to full emissivity (re = 1) 
agrees very well with data within the experimental error bars (shown in the plot), whereas the 50% 
emissivity case is completely excluded by data, since the ratio is too high in that case, reflecting 
what we saw in the two previous plots. The values of the parameters are the same as before, i.e., 
T = 411 MeV, T c = 175 MeV, T f = 150 MeV, and the transverse radius, R T ~ 7 fm/c. 
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